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We describe how to compute planar gluon scattering amplitudes at strong coupling in 
Af = 4 super Yang Mills by using the gauge/string duality. The computation boils down 
to finding a certain classical string configuration whose boundary conditions are determined 
by the gluon momenta. The results are infrared divergent. We introduce the gravity version 
of dimensional regularization to define finite quantities. The leading and subleading IR 
divergencies are characterized by two functions of the coupling that we compute at strong 
coupling. We compute also the full finite form for the four point amplitude and we find 
agreement with a recent ansatz by Bern, Dixon and Smirnov. 



1. Introduction 



In this article we describe a method for computing gluon scattering amplitudes at 
strong coupling in A/" = 4 super Yang Mills. Of course, these amplitudes are infrared 
divergent and are not good observables. Nevertheless, in practical computations for collider 
physics it is often useful to compute these amplitudes as an intermediate step towards 
computing actual well defined observables. Proper observables are IR finite (IR safe), see 
e.g Furthermore, in QCD computations it has proven useful to know the A/" = 4 

super Yang Mills result as a building block [Q. Thus, it is interesting to understand the 
behavior of scattering amplitudes at strong coupling. We perform the strong coupling 
computation by using the gauge theory /gravity duality that relates Af = 4 super Yang 
Mills to string theory on AdS^ x We consider planar amplitudes. On the string 

theory side, the leading order result at strong coupling is given by a single classical string 
configuration associated to the scattering process. A similar result was found by Gross 
and Mende in their investigation of fixed angle, high energy scattering of strings in fiat 
space. As we explain below, the string theory scattering in AdS is happening at fixed 
angles and large energy and it is thus determined by a classical solution. The final form 
for the color ordered planar scattering amplitude of n gluons at strong coupling is of the 
form 

where Sd denotes the classical action of a classical solution of the string worldsheet equa- 
tions, which is proportional to the area of the string world-sheet. The solution depends on 
the momenta, /c^, of the gluons. The whole dependence of the coupling is in the overall 
factor. Of course, we expect corrections which we do not compute. The structure of 

the IR divergences is precisely as expected in the field theory. This comparison enables us 
to compute the strong coupling expression for the function QoW EJl'SHMl characterizing 



the subleading IR divergent terms. The IR regularization is done via dimensional regular- 
ization, as in the field theory. This is achieved by using the gauge theory/gravity duality 
for Dp-branes for general p and then performing an analytic continuation in p = 3 — 2e. 

One of the motivations of this work was the very interesting conjecture by Bern, 
Dixon and Smirnov (see also |T^) for the all order form of the n gluon MHV scattering 
amplitudes. We have computed explicitly the full strong coupling answer for the four point 
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amplitude and found precise agreement with their conjecture. Their conjecture is that the 
four point amphtude has the form 



exp 



fix) 

(IR divergent) H — (fog(s/t))^ + (constant) 
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where s, t are Mandelstam variables, /(A) is directly related with the cusp anomalous 
dimension and the IR divergent terms are well characterized Sudakov-like factors 0. 

Scattering amplitudes via the AdS/CFT duality have been study in many articles. 
See for example [12, p!3|jn ,15,16,rf| and references therein. 

This paper is organized as follows. In section 2 we motivate and describe the general 
prescription for computing the leading order approximation for the amplitudes. In section 
3 we compute explicitly the classical string solution that describes the four point amplitude, 
then we proceed with a discussion of the structure of infrared divergencies and perform the 
comparison with the results in 0. Some conclusions and open problems are sketched in 
section 4. Many additional remarks and technical details are included in several appendices. 



2. Gluon scattering amplitudes 

Since scattering amplitudes of colored objects are not well defined in the conformal 
theory it is necessary to introduce an infrared regulator. The answer we obtain will de- 
pend on the regularization scheme. Once we compute a well defined (IR safe) physical 
observable the IR regulator will drop out. An example of a well defined physical observ- 
able is the amplitude for a process involving narrow jets going in some specific angular 
directions. The answer will depend on the precise definition of the jet observable, but not 
on the IR regulator. Further discussion of these issues can be found in A popular 

regularization scheme is dimensional regularization and we will use it in the next section. 
For the time being, it is convenient to use a different IR regulator. In terms of the gravity 
dual, this IR regulator is a D-brane that extends along the worldvolume directions but is 
localized in the radial direction. In other words, we start with the AdS^ metric 



(2.1) 



^2 



and we place a D-brane at a large value zifi. In terms of the field theory, such D-branes 
arise, for example, if we go to the Coulomb branch of the theory. The asymptotic states 
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are open strings that end on the D-brane0. We then scatter these open strings. We are 
interested in keeping the momentum fixed as we take away the IR cutoff. This means 
that the proper momentum of the strings hving at zj^ is very large, kzi^ ^ 1, where k 
is the momentum in field theory units {k is conjugate to translations in x). Thus, we are 
studying the scattering of open strings at fixed angles and very high momentum. Such 
amplitudes were studied in fiat space by Gross and Mende, 0. The important feature 
noted in [0, is that amplitudes at high momentum transfer are dominated by a saddle 
point of the classical actioni. Thus, in order to compute the amplitude we simply have 
to compute a solution of the classical action. In our case we need to consider a classical 
string in AdS. 




s 

Fig. 1: Order of the particles in the diagram and definition of s and t for the 
four point amplitude. 



Let us first recall the fiat space case [§. We consider scattering of open strings on a 
D brane that is transverse to the direction z, sitting aX z = zm. We want to consider a 
worldsheet with the topology of disk with vertex operator insertions on its boundary. See 
figure fig. 1. Near each vertex operator insertion the solution should behave as i 

ik^" log \w\'^ (2.2) 

^ The color structure of the amplitude is not very apparent with this regularization. We can 
get a ghmpse of it by introducing A'' branes in the IR (instead of just one). In other words, we 
go to a generic point in the Coulomb branch. Then the external gluons are strings stretching 
between these branes. The color ordered amphtude would arise when we consider a configuration 
where one open string ends on the brane where the next open string starts. Namely, the ith open 
string goes between branes i and i + 1. 

^ Here we concentrate on genus zero amplitude, Q computed also higher genus amplitudes 
using the same idea. 

^ We set a' = 1. 
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where k'^ is the momentum of the open string. In addition we require z = zm on the 
boundaries of the worldsheet. The solution is 



= 1 



i'^k'^\og\w - Wi\^ , z = ziR (2.3) 



ActuaUy, we stiU need to determine Wi. The actual solution that minimizes the action has 
fixed values of Wi, up to conformal transformations on the worldsheet. For example, in 
the case of the four point function we can set three of the Wi to wi = 0, ws = 1, 104 = 00 
and then W2 is determined by inserting the above solution into the action and minimizing 
with respect to W2 which gives 

W2 = ^-, s = -{ki+k2)\ t = -{k^+k4f (2.4) 

s + t 

The value of the action is then 

Sfiat{s, t) = s log(-s) + t log(-t) -{s + t) log(-s - t) (2.5) 

and the leading approximation to the scattering amplitude is then A = . Notice that 
the worldsheet has Euclidean signature in the regime that s, t < (spacelike s and t 
channel momentum transfer, and timelike u channel momentum transfer). In this regime 
the solution ( |2.3| ) is such that the coordinates along the brane, x^, are purely imaginary. 
This amplitude is exponentially suppressed and it represents the very small probability 
process where the two incoming string states tunnel to the two outgoing ones. Notice 
that we have an open string amplitude and the order of the boundary vertex operators 
is important. Note that this leading exponential behavior (|2.5| ) is independent of the 
particular string states we are scattering, as long as we keep them fixed when we increase 
the momentum transfer. 

Let us now consider the AdS case. Now it will be more difficult to find the classical 
solution. However, there is one important aspect of the classical solution that we can 
understand in a qualitative way. Namely, we expect that the solution will be such that 
in the central region of the collision the value of z will be much smaller than zir and 
that it will be roughly proportional to the inverse of or since they set the off 
shell momentum transfer of the process. This expectation is based on the idea that the z 



direction should correspond to the off shell energy scale of the process. In |14] a similar 



scattering process was considered and it was also found that the leading contribution came 



from 1/z ~ a/ (—■§)• ^ 

In order to state most simply the boundary conditions for the worldsheet it is conve- 
nient to describe the solution in terms of T-dual coordinates defined in the following 
way. We start with a metric that contains 

ds^ = w^{z)dXf,dxf' + ■■■ (2.6) 

where w is the warp factor. We define T-dual variables by 

day'' = iw''iz)e^pdpx'' (2.7) 

In the regime under consideration the T-dual coordinates are real and the worldsheet is 
Euclidean. In addition, the boundary condition for the original coordinates x'^, which is 
that they carry momentum k'^, translates into the condition that y^ has "winding" 

Ay" = 2TTk" (2.8) 

Note that we are not taking the coordinates to be compact (specially time!). One can 
view this as purely a mathematical operation that makes it easier to find the classical 
solutions.i The T-duality will also produce a non-trivial dilaton field, but it will not affect 
our classical solutions. The T-dual metric is again AdS^ after defining r = — i 



ds^ = R 



2 



dy^dy" + dr^ 
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r = — (2.9) 

z 



The full Green Schwartz string action was computed in these "T-dual" variables in |]T^ . 



The main difference with our configuration is that in |14] the asymptotic states were closed 
strings in the bulk. Thus they could move in the z direction more easily than the open strings we 
consider, which are attached to a D-brane at zir. The open strings can move in z by stretching 
away from the D-brane, which is a stringy excitation. 

^ One could consider situations with compact spatial coordinates. Then we would be doing 
an honest T-duality on the spatial coordinates. 
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Notice that we do not do a T-duality in the z direction. 
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Fig. 2: The kinematic data, namely the sequence of momenta k'^, - ■ ■ ,kn, trans- 
lates into a sequence of lightlike segments joining points in the T-dual space 
parametrized by y'^. These points are separated by 27rk^ . This curve lives at 
r = in the T-dual AdS space ( |2.9| ). We need to find a minimal surface that ends 
on this curve. 



The solution we want is a surface which at r = R'^/zin ends on a particular one 
dimensional line which is constructed as follows. For each particle we have a lightlike 



segment joining two points separated by (|2.8|) . We concatenate these segments according 



to the ordering of the open strings in the disk diagram. This ordering is interpreted as the 
particular color ordering of the amplitude. See fig. 2. The resulting line consists of lightlike 
segments. The condition that the line closes corresponds to the momentum conservation 
condition. As we explained above the solution lives at values of r > rjfi = R'^/zm. As 
we take the limit ziji — oo we find that the boundary of the worldsheet moves to the 
boundary of the T-dual metric ( pj.9| ) which is at r = 0. From the point of view of the 
T-dual metric ( ^.9| ) the computation that we are doing is formally the same as the one we 
would do [|T^ if we were computing (in the classical string approximation) the expectation 
value of a Wilson loop given by a sequence of lightlike segment J3. 

In conclusion, the leading exponential behavior of the n-point scattering amplitude is 
given by the area. A, of the minimal surface that ends on a sequence of lightlike segments 
on the boundary of (|2.9|) 

^^e-^^ = e-^^ (2.10) 



The area A = A(/c^, • ■ ■ , /c^) contains the kinematic information about the momenta. The 
amplitude is color ordered in a way that is reflected by the particular order in which the 
lightlike segments are arranged along the boundary. There is no information about the 
particular polarizations or states of the gluons, which contribute to prefactors in ( |2.10|) , 
and are of subleading order in 1/ y/X. From the worldsheet point of view, we will need to 



^ This prescription is vaguely reminiscent to the one used for non-commutative theories 
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supply this information once we quantize around the classical solution. Such terms are 
beyond the scope of this paper. Note that the A dependence is contained purely in the 
factor multiplying the areai. 

The result ( |2.1C1| ) is still somewhat formal since the area is infinite due to the infrared 
divergences that we mentioned above. In order to find a finite answer we will need to 
regularize the result. In the next section we will discuss the structure of the infrared 
divergences and we will do an explicit computation for the four point amplitude. 

A reader familiar with AdS / CFT might be surprised by the fact that in terms of the 
original coordinates we are setting boundary conditions at r = or 2; = oo , rather than 
at 2; = 0. This apparent confusion goes away once we consider AdS in global coordinates, 
which is better in order to see whether we are at the boundary or not. In those coordinates 
one can see that the surface z = 00 indeed intersects the boundary, as a;^ — 00 and we will 
later check explicitly that the solution in terms of the coordinates {x^,z) is intersecting 
the boundary of AdS. This is discussed in more detail in appendix A. Thus, there is no 
contradiction with the general principle stating that good observables are defined on the 
boundary of AdS. 



3. Computation of the four point amplitude at strong coupling 

In this section we consider the planar four gluon amplitude. We label the momenta 
as /ci, /c2, ks, /c4, where the subindex indicates the color ordering, see fig. 1. We consider 
the region where particles 1 and 3 are incoming and particles 2 and 4 are outgoing. We 
label by k the center of mass energy or momentum of each of the incoming particles and 
we denote by <f> the scattering angle in the center of mass frame. We introduce the usual 
Mandelstam variables 

s = -{ki + k2f = -2ki.k2 = -Ak^ sin^ | 

t = -{ki + k^f = -2ki.ki = -Ak^ cos^ I (3.1) 
u = -{ki + k^f = -2ki.ks = Ak^ = -{s + t) 

We will focus on the region where s,t < (they correspond to spacelike momentum 
transfer) . 

^ The area is computed in an AdS space with radius one. 



7 



Fig. 3: Sequence of lightlike segments which specifies the scattering configuration. 
This figure lives at r = of the metric ( p.QI ). 



As we explained above, we should find a classical string solution specified by these 
momenta. It is simplest to think about the solution in T-dual coordinates where the 
problem boils down to finding a minimal surface in the T-dual AdS^ space ( |2.9| ). This 
surface ends at r = on a closed sequence of lightlike segments whose sides are specified 
by the lightlike momentum vectors (27r)/cf , see fig. 3. 

3.1. The lightlike cusp 

We start by considering the solution near the cusp where two of the lightlike lines 
meet. So we consider two semi infinite lightlike lines meeting at a point. This case was 



considered in |21] and it will prove useful for generating the solution we want. It is a 



surface that can be embedded in an AdSs subspace of AdS^ 




Fig. 4: The lightlike cusp. The thin lines shows the light cone. 



8 



We are interested in computing the surface ending on a light-like Wilson loop which is 
along = ±y^, > 0,i see fig. 4. The problem has a boost and scaling symmetry that 
becomes explicit if we choose the following parametrization 

yo = cosher, yi = sinhcr, r = e^w (3-3) 

Boosts and scaling transformations are simply shifts of a and r. Then the Nambu-Goto 
action becomes 



One can explicitly check that w{t) = \pi solves the equations of motion, hence the surface 
is given by 



v^V 2/0 - 2/? = V2 V^rt^ , y±=l/0±yi (3.5) 



When we insert the solution in the action ( pi. 41 ) we find that the lagrangian is purely 
imaginary, this means that the amplitude A ~ e*"^ will have an exponential suppression 
factor 

iS = —Se = — / dadr (3.6) 



An _ 

where 5" is the action for a spacelike surface embedded in the Lorentzian target space that 
we are considering. This integral is infinite. We will later discuss its regularization. 

It is instructive to study this solution in terms of embedding coordinates. These are 
coordinates where we view AdS^ as the following surface embedded in R^'^ 

-y^i - + y2 ^ Yi + Yi + Yi = -i (3.7) 



The relation between these and the Poincare coordinates in (|2.9| ) is 



Y^^!^ , ^ = 0,---,3 
r 



(3.S 



^ One can also consider Wilson loops along y° = =by^, > 0, the basic difference with the 
ones considered here is that their world-sheet is Lorentzian and z is imaginary. 
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We can now write the surface corresponding to the cusp in terms of the equations 



21 



Y^' - = Yt - Yl, ^2 = I3 = (3.9) 

3.2. The four lighlike segments solution 

We now consider a Wilson loop containing four light-like edges, which contains four 
cusps like the one considered above. The configuration of lightlike lines is shown in fig. 3. 

In order to write the Nambu-Goto action it is convenient to consider Poincare coor- 
dinates (r, yo, yi, 2/2), setting 7/3 = 0, and parametrize the surface by its projection to the 
(2/1,2/2) plane. We consider first the case with s = t where the projection of the Wilson 
lines in fig. 3 is a square. The Nambu-Goto action is then the action for two fields yo and 
r living on a square parametrized by yi and 2/2- The action reads 



e [ J J v^l + (a,r)2 - {diyoY - (91^^22/0 - d2rdiyo)'^ 

S = — dyidy2 5 (3.10) 

An / 



By scale invariance, we can change the size of the square. We choose the edges of the 
square to be at 2/1,2/2 = ±1- The boundary conditions are then given by 

r(±l, 2/2) = r(2/i, ±1) = 0, 2/o(±l, 2/2) = ±2/2, 2/0(2/1, ±1) = ±2/i (3.11) 

Prom the solution for the single cusp, we can obtain, after boost transformations, the form 
of the solution in the vicinity of any of the cusps. The following expression for the fields 
can be easily seen to have the right behavior close to the cusps 



2/0(2/1,2/2) = 2/12/2, r(?/i,2/2) = y(l-2/?)(l-2/i) (3-12) 

Remarkably it turns out to be a solution of the equations of motion. When expressed in 
terms of embedding coordinates, the surface is given by the equations 

Ys^O, n^O, YoY_^=YiY2 (3.13) 

In fact this solution is related by the AdS isometries {SO{2, 4) transformations) to the cusp 
solution ( The reader might be puzzled by the following. We seem to have mapped a 

solution with two lightlike lines on the boundary to one with four lightlike lines. The AdS 



10 In order to go from to ( p^ ) we take ¥2 ^ Y4, Yo ^ -^{Yo+Y^i), F-i ^ 75 (^o-^i^-i), 
^1 ^ 75(^1 + ^2) and y4 ^^{Yi- Y2) 
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isometries are conformal transformations on the boundary and conformal transformations 
preserve angles and cannot produce cusps where there were none. The solution to this 
apparent puzzle is that the cusp solution ( ^^91) really has four cusps once it is embedded 



in global coordinates pTJ. We miss some of the cusps when we use Poincare coordinates 
because those cusps are on the boundary of the Minkowski space parametrized by y^. In 
fact this is a simpler alternative way to derive the solution. Namely we start with the 
cusp solution (|3.9|) , notice that it really has four cusps and then map it through conformal 
transformations to the solution we really want ( |3.12| )-( ^.13|) . 

After we understood this point it becomes a simple exercise to compute the solution 
for general s and t. We simply need to apply an 5*0(2,4) transformation to the solution 



we already have. Starting from (|3.13|) we perform a boost in the 04 plane and obtain 



Y^-vYo = Q, y_i7(yo - t^n) = 7"'>oy-i = , l3 = (3.14) 



where 7""*^ = yl — v^. 

Let us now write the solutions in terms of worldsheet coordinates in conformal gauge. 
Let us first go back to the solution for the case with s = t, ( |3.12| ), and compute the induced 
metric on the worldsheet. We find 

^•^^ —7 tTo' + 7 ^9^0" = ^"^1 + du\ , where yi — tanhu^ (3.15) 

[l-ytY (I-2/2) 

Notice that the metric on the worldsheet is Euclidean. More precisely, we have a spacelike 
surface embedded in a Lorentzian target space. Written in terms of Ui coordinates the 
solution ( |3.12| ) becomes 

yi = tanhwi , y2 = tanh'U2 , r = \ \ , yo = tanhwi tanh'U2 

cosh Ml coshM2 

Yq = sinh-ui sinhw2 , ^1 = sinhwi coshu2 , Y2 — coshui sinhu2 , 

y_i = coshwi coshM2 , I4 = Y3 = 

(3.16) 

This is now a solution of the equations in conformal gauge, whose action reads 

(g-^ + W (3.17) 

2n J 2n J 2 
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Note that the metric ( 3.15| ) is Euchdean, this is responsible for the extra i in this formula. 
The lagrangian density evaluated on the solution is simply C = 1. Performing the boost 
( p.l4|) and a simple rescaling we now find the solution for s ^ t 



a 



yi 



cosh ui cosh U2 + b sinh ui sinh U2 

a sinh Ui cosh U2 
cosh ui cosh U2 + b sinh ui sinh U2 



yo 



2/2 



aVl + b"^ sinh Ml sinh 



U2 



cosh ui cosh U2 +b sinh ui sinh U2 
a cosh Ml sinhu2 



cosh ui cosh U2 +b sinh ui sinh U2 

(3.18) 

where b = and we consider 6 < 1. The parameter a sets the overall scale of the 
momentum. The solution approaches the boundary of AdS^ where ui or U2 go to plus 
or minus infinity. These four possibilities correspond to the four lightlike lines on the 
boundary. For example, if we take ui +00 we find that r = and 

a atanhw2 / 7:7 tanhw2 



yi 



l + 6tanhM2 ' l + 6tanh'U2' v ■ - l-|-5tanhu2 

We see that yi + by2 = a and that we have a lightlike line going between two points whose 
projections on the 1/1,2/2 plane are located at 



yo 



iVl + b^- 



(3.19) 



A : yi=y2 
which are reached at U2 ±00. 



1 + 6 



and B : yi = -y2 



a 



1-6 



(3.20) 





(b) 

Fig. 5: Projection of the light like lines on the yi,y2 plane for (a) s = t and (b) 
s ^ t. The line also moves in the time direction with a slope such that we get 
a lightlike line. Points on opposite vertices sit at equal times. Time goes up and 
down as we move from segment to segment. 



By considering other limits we get the other segments. The values of s and t are given 
by the square of the distance between the vertex of two non-adjacent cusps, see fig. 5. In 
terms of the parameters a and b they are given by 
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where the factors of 27c comes from 

The solution and the value of the action are symmetric under s ^ t, which is a 
symmetry of the full problem. 

It is also instructive to write the solution in terms of the original AdS coordinates 
(^). We obtain 

iR^ /ii2 1 , /""ill 
xi = A sinh2w2 + b{-— + - sinh2wi) 

X2 = I -— smh2ui + b{— smh2u2) 



iR 



2 



(3.22) 



xq = V 1 + 6^ (cosh M2 — cosh Ui) , 

2a 

z = — (cosh Ml cosh'U2 + bsinhui sinh'U2) , 
a 

We see that for large Ui the solution is such that it carries momentum in the spatial direc- 
tions. The solution lives in complexified AdS space since it represents a sort of tunnelling 
solution. For Ui ±oo the solution goes to the region where z is large, which naively seems 
to corresponds to the IR of the field theory. On the other hand, we are also finding that x^ 
are going to infinity at the same time. We can go to global coordinates and find that the 
solution ( p.22|) indeed approaches the boundary of AdS. It touches the boundary on the 
surface with 1/^ = which corresponds to the region of the boundary that is the boundary 
of the Penrose diagram of the Minkowski slices parametrized by x^. See appendix A. In 
the central scattering region, Ui ~ 0, we have that z ~ R^ /a ~ -R^("y= + :^7=t-' "^^^i^h is in 
agreement with the discussion in section two, where we said that the scale of momentum 
transfer should set the minimum value of z that the solution explores 0. 

We are now almost ready to evaluate the action. However, we should note that there 
is a small subtlety. The action in terms of the original coordinates and the action in 
terms of the T-dual coordinates differ by a total derivative which will contribute with a 
boundary term. The correct action to evaluate turns out to be the area of the surface 
in the T-dual coordinates. When we consider the problem in the original coordinates we 
should remember that we are putting boundary conditions that fix the momentum at each 



The above behavior can be quahtatively understood as the result of two competing factors. 
On one side, high proper momentum transfers suppress the scattering and this suppression can be 
lowered by going to regions of smaller z. On the other hand we have to pay a price for stretching 
the string away from z = zir in the z direction. 



13 



boundary. In order to have a proper variational principle we should add a boundary term. 
This boundary term is precisely the one that turns the original action into the T-dual 
action. 



3.3. Dimensional regularization in the gravity dual 

As we mentioned above the action of the solutions we have discussed is infinite. A 
popular regularization method for A/" = 4 super Yang Mills is the so called dimensional 
reduction scheme |^^. In this scheme one goes to a general dimension D = 4 — 2e but 
one continues to use a theory with 16 supercharges. In other words, one considers the 
dimensional reduction of ten dimensional super Yang Mills to 4 — 2e dimensions. For 
integer dimensions these are precisely the low energy theories living on Dp branes, where 
p — D — 1. The gravity dual of these theories involves the string frame metric 

ds^ =r^'^dxl + f^'^[dr^ + r^dnl_j,] , D = 4 - 2e 

, CpXp r,4e^3e-r('o , \ (3.23) 

/=7CT ' CD =2 TT r(2 + e) 



where Xd = gf^N and g^ is the couplinglB. We parametrize the coupling in D dimensions 
in terms of the IR cutoff scale fi as it was done in the field theory analysis of 

= ^ = (^-^^^ 

where A = A4 is the dimensionless four dimensional coupling which is kept fixed as we vary 
e. 

From now on, we will drop the sphere part of the metric since it does not play an 
important role. We now compute the metric for the T-dual variables defined through ( |2.7| ) 
(with w'^ = /"^/2). We get 

ds' = f'l\dyl + dr^) = (^^l^^) (3.25) 

Notice that the IR region of the original metric corresponds to the region where r ~ in 
the T-dual metric {^^. 

In the region very close to r ~ we cannot trust the gravity dual and we should use 
the weakly coupled field theory in the IR. In this field theory picture we can define the 
asymptotic gluon states. We will see that we will only need the geometry in the region 
where we can trust it if we are interested in the strong coupling behavior. 



This is defined as the coefficient of the D dimensional action in the usual way S = -^Tr[F'^] 



and rr[r"r''] = Uab where are the generators of SU(N). 
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3.4- Evaluation of the action in dimensional regularization 

We regularize the theory by considering L>p-branes, with p = 3 — 2e. We are then lead 
to the following action 



S = [ ^ (3.26) 



where £e=o is the lagrangian density for AdS^, as in ( p.lO| ) or ( |3.17| ) 



In order to understand how dimensional regularization works, let us perform the 
computation for the lightlike cusp in general dimensions. We still have the boost symmetry 
and we can make an ansatz similar to the one in (^^) but now the Lagrangian depends 
explicitly on r 

S = [da[ dre- ^^'^^J^'^ (3.27) 

It turns out that a constant w is still a solution, but now the constant is 



w = \/2^1 + ^ r = v^Vl + e/2v/2/+2/- (3.28) 

Inserting this solution into the action and writing the integral as an integral over we 
get 



{2y+y-y+^/' e^2y+y^y/^' 87r(l + e/2) 1+^/2 ^^'"^^ 

where y^ is a cutoff for large y^ 0. We will later see that this double pole is in agreement 
with field theory expectations. One might be worried that the metric (|3.25| ) is getting 
highly curved as r — > (for e < 0). Fortunately this is not a problem if we are only 
interested in determining the strong coupling behavior of the amplitude. We see this as 
follows. The integral is cutoff at the point where 

~ 1 (3.30) 

But the effective curvature, TZ, at that point is of the form 

n ~ < 1 , if A > 1 (3.31) 

V A 
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In the computation of the actual amplitudes this will be the momenta of the particles. 
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Thus, for large lambda we can perform dimensional regularization without worrying about 
the region with strong curvature. If we wanted to understand the result at all values of A, 
then it would be important to understand the whole region. A simple way to put it, is to 
say that dimensional regularization and strong coupling commute. 

Let us now turn to the problem of the four-point scattering amplitude. When e 7^ 
we have a different lagrangian ( p.26 ) and we would need to find the solutions for the new 



lagrangian, as we did above for the cusp. Fortunately, there is a simple trick that allows 
us to find the solution to the accuracy that we nee d0. We first note that if we have a 
lagrangian which has the expansion L = Co + eCi +€^£2, then we can expand the solutions 
of the equations of motion as g = go + e?i + ■ ■ ■• We will be interested in evaluating the 
final answer to zeroth order in e. However, since we have IR divergencies we find that 
the leading order solution qo will give a leading double pole in e. Thus when we do the 
formal expansion we mentioned above, we will want to evaluate the action to the formal 
order e^, which in reality will be order e''. Note that then we will not need to know the 
solution to second order in e since that will contribute a term of the form J | q2 which 
vanishes due to the fact that qq obeys the zeroth order equations of motion. For a similar 
reason the first order solution qi will only contribute to terms of formal order e^. Thus 
to real orders and we can simply evaluate the zeroth order solution in the new, 
e dependent lagrangian and we will get an accurate enough answer. Since the first order 
solution can only contribute to formal order we will need an IR divergence of order l/e'^ 
to give a finite answer. These divergences only arise in the cusps. Thus, we only need 
the gi solution near the cusps. However, for the cusp region we know the solution ( |3.28|) . 
When we have the general zeroth order solution of the problem, we can get a solution that 
is accurate enough at the cusps by writing 



1 + I r,=o , ~ y^=o (3-32) 

where the e = solutions are the ones we discussed above. Inserting these expressions 
into the action ( p.26|) gives us an accurate enough answer to extract the finite pieces in 
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In principle, one could compute the scattering amplitudes for gluons in the case of D = 5, 6 



if the center of mass energy is such that we can trust the corresponding gravity solutions [23|. 
That would require solving the equations for the new lagrangian (|3l26|) with e = 
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the amplitude c3. A more detailed analysis of the finite e equations in the various regions 
(near the cusps, near the lines) shows that the above argument is indeed correct. 
By inserting these expression into the action we get 

/oo 
duidu2{cosh.ui cosh W2 + 6sinhui sinhw2)'^ (l + e/i + e^(/2 — 2/|) + ...) 
-oo 

(3.33) 



where 



^ _(62 - l)(cosh2ui + cosh2u2) - 2(1 + 6^) 
^ 8(coshwi cosh W2 + fesinhwi sinh ^2)^ 

^ 1 + 6^ — (1 + 6^) cosh2ui cosh2w2 — 26sinh2'Ui sinh2w2 (3 34) 

^ 16(coshui coshM2 + 6sinhui sinhu2)^ 

„ VXdCd 1 
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where we have expanded up to terms that give finite order answers in the final result. The 
integrals can be performed as explained in appendix B. The final result is 

= -5. 2F,(^, b^) + 1/2) (3.35) 

It is then straightforward to expand in powers of e up to finite contributions. We need to 
recall the expressions for s,t in ( 3.21 ) and also the formulas for cd ( |3.23| ) and Ad ( |3.24| ). 
Putting all this together we get the final answer 



e'^ ^ exp 



iSd^v + ^ (log -Y + C 



(3.36) 



where Sdiv,s and Sdiv,t are the divergent pieces associated to each cusp or pair of consecutive 
gluons. There are two pairs with total squared momentum t and two with s. We have 



and Sdiv,t is given by a similar expression with s t. We now compare what we obtained 
here with the field theory results and conjectures in 0. 



In order to compute the regularized area we use the Nambu-Goto action in the static gauge 
chosen above, since the corrected solution near the cusps ( ^.321 ) is expressed most simply in this 
gauge. 

We have started with a cutoff scafe as done in the field theory analisys of 0, however, one 
can verify that if we have started with y^/x instead, all the terms containing log 2 in our final 
result would disappear. 
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3. 5. IR divergences of amplitudes 

In this section we recall some general results on IR divergences of scattering amplitudes 
0,11,^. We will focus here on planar diagrams and color ordered amplitudes. The first 
result is that the IR singularities of the amplitude can be associated to consecutive gluons in 
the color ordered amplitude. This is fairly clear once we recognize that the IR singularities 
come from low momentum gluons and we use that we consider only planar diagrams. The 
leading divergence goes like 

^ = e-^(iogM)^ (3.39) 

where /U is a mass scale that is acting as an IR cutoff. We have one factor like ( |3.3y| ) for each 
pair of consecutive gluons. The reason that the divergencies exponentiate in this way is the 
following. The divergencies come from the exchange of soft gluons among two consecutive 
hard gluons. In the limit that the hard gluon momenta are infinite, the configuration of 
hard gluons is invariant under two symmetries: boosts and scale transformations. Let us 
denote by a the boost parameter and by r the parameter generating scale transformations, 
x'^ e'^x^. Thus naively the amplitude will be of the form 

Ad^. = e-'^(^)^-^- (3.40) 

where Act and Ar is the range of scale and conformal transformations where the approx- 
imation of exact boost and scaling symmetry are valid. (Recall the classical string result 
( p.6|) ). This is conceptually similar to the problem of computing the partition function of 
a system that is invariant under time translations and spatial translations. The answer 
will he Z = e~^^'^ where L and T are two IR cutoffs. In fact, this can be made more 
explicit by choosing coordinates in R^^^ that lead to a metric which is Weyl equivalent 
to a metric where these symmetries are explicit^. The function / that appears in ( |3.3y| ) 
is proportional to the cusp anomalous dimension for a Wilson loop in the fundamental 



representation li^. Of course we also need another scale so that the log in (|3.39|) makes 



One can write the metric of R^^^ as ds^ = e^'^{—dT^ + ds^^). Since the theory is conformal 
invariant we can drop the conformal factor in the metric. We can further choose coordinates in 
hyperbolic space so that we have now the metric ds = —dr + dp -\- cosh pda + sinh pd(p . 
Boosts and scale transformations coorrespond now to shifts of r and a. 
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For two Wilson lines forming a spacelike cusp where the two lines differ by a large boost 
parameter 7 the anomalous dimension is {W) ~ ^-^~i^°s.i^iR/ ^uv) ^ where L are UV and IR 



cutoffs. We refer the reader to [^,24,25] and references therein for a discussion of all these ideas. 
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sense. This is provided by —s = {pi +^2)^ = 2piP2 where pi and p2 are the momenta 
of the two hnes. The function / in (|3.39| ) also appears when one computes the dimension 
of operators of high spin, S, of the schematic form Tr[$9'^$], with $ in the adjoint |25 . 
These operators have twist 

A-S = fiX)\ogS (3.41) 



This function was computed perturbatively up to four loops in ||2^,^ , at strong coupling 
using strings in AdS in [^|2^ and given exactly as a solution of an integral equation in 
Il30[] using integrability (see also [0). Of course, in addition to the leading divergence in 
( p.39|) we can have a subleading divergence involving a single log. Thus, we can introduce 
a second function g through 

Aa..,s = exp ^ -4^ f log -j^) ' - ^ ( log ) ^ (3.42) 




16 V (-s). 

Where we have defined a new function g. The precise form of g will depend on the precise 
definition of the IR regulator since shifting the log in the first term by a constant can affect 
the form of g. In other words, changing fj, ^ we change 

g^g + flogK (3.43) 

Let us now review how these divergences appear when we perform dimensional reg- 
ularization. The double logs in ( p.39|) will arise if we have a double pole in e. Thus, in 
dimensional regularization the divergences should organize as 

.4.„,.^e.p|-i./<-.f^Vi.'-'>mi (3.44) 



This formula, together with (|3.24|) , gives a precise definition for g. We now see that 
expanding in e we reproduce the double logs in (|3.39| ) if 

(a^) /^-^)(A) = /(A), X^g(-'\X) = g{X) (3.45) 

By comparing the general expression for the IR divergence ( |3.44| ) to our result ( |3.38D we 
can compute the functions / and g at strong coupling. We find 

/ = 4^, ,= ^2(l-log2) (3.46) 

Notice that the square root of A dependence introduces a factor of 2 when we integrate as 
in ( |3.45|) . Of course, the function / has the same value that was computed in other ways 
in [p8| , [2T[] this is implied by the general theory of IR divergences we have just reviewed. In 
appendix C we discuss the extrapolation of the weak coupling results to strong coupling 
and the comparison with our answer ( |3.46|) . 



What we called g here is called Qo in 
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3. 6. Field theory results for the four point amplitude 

The four point scattering amplitude in maximally supersymmetric theories in D di- 
mensions has the form 

A = ATreea{s,t) (3.47) 

This is completely determined by the fact that the theory has 16 supersymmetries. S A 
simple way to understand it is the following. In a theory with 16 supersymmetries the 
generic massive representation has 2^ components. The gluons live in a representation 
which has 16 = 2^ components. This is possible because it is a massless representation 
and half of the supercharges act trivially. When we think about a 2 — > 2 scattering process 
we can see that the state formed by two gluons transforms as a massive representation and 
has 2^ states, which is the same as the total number of polarization states for two gluons. 
This means that there is a unique intermediate state. This fixes the S matrix uniquely 
up to a common function.!! Of course this S matrix determined by the symmetries will 
be equal to the tree level S-matrix which preserves all symmetries. This explains the form 
of the four point scattering amplitude in ( |3.471 ) in all dimensions. In four dimensions the 
function a{s,t) has to obey further constraints. If we demand scaling symmetry we would 
conclude that it is a fimction of a{s/t) and if we further demand that it is invariant under 
special conformal symmetries we would conclude that a is a constant. However, due to 
the IR divergencies we have an additional dependence on the cutoff, which enables the 
scattering amplitude to be a non-constant function of s, t. HI 

Bern, Dixon and Smirnov made a conjecture for the exact form of the four loop 
amplitude of the form 

A = Atree{Ad^.,sf{Ad^.,tf exp |^[(log jf + AtI^ /?,] + C(A)| (3.48) 

where C(A) is only a function of the coupling, / is the same function appearing above ( |3.39|) , 
and Adiv,s is given by (|3.44|) . We see that the momentum dependent finite piece of our 



strong coupling expression has precisely the form predicted by [0 , including all numerical 
factors after we include the appropriate strong coupling form for the cusp anomalous 



See for instance [32,33,34] and appendix E of |35]. 
A similar argument was used in [36| for S matrices on spin chains. 

It seems possible that the functional form of the IR regularized result is still determined by 
the action of the special conformal generators, though we were not able to prove it. 
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dimension (|3.46| ). Unfortunately we cannot test 0] for the the constant finite pieces, 
though we could do it if we computed the n point amplitudes, which are also predicted in 
. Of course, what we can definitely say is that ( |3.37| ) implies a strong coupling value for 
a combination of the constants introduced in 0. 



4. Conclusions 

In this paper we have given a prescription for computing planar gluon scattering 
amplitudes at strong coupling using the gauge theory /gravity duality. The computation 
involves finding a classical string solution moving in (complexified) AdS^ with a prescribed 
asymptotic behavior determined by the gluon momenta. The computation can be regu- 
larized using the gravity version of dimensional regularization. Though our prescription 
would work for an arbitrary n point amplitude, we have only computed explicit answers 
for the four point amplitude. We found results that agree in all detail with the conjecture 
of Bern, Dixon and Smirnov . The structure of IR divergences is precisely as expected 
from general field theory reasoning. The detailed comparison enables us to extract the 
strong coupling form of the function that determines the subleading divergent terms. The 
leading divergent terms are determined by the cusp anomalous dimension which is already 



known at strong coupling |2S][21 



It is amusing that the computation is mathematically similar to the computation we 
would do in order to compute the expectation values of locally lightlike Wilson loops. This 
formal similarity, however, might disappear at higher orders in 

It seems possible to try to check the conjecture of [0 for n point functions for n > 4. 
Since the AdS^ sigma model is integrable, it should be possible to find the appropriate 
classical solutions and evaluate the action. Or even evaluate the action without finding 
the full explicit form of the classical solutions! 

It would also be interesting to learn how to do computations of scattering amplitudes 
exactly as a function of the coupling by using integrability and spin chains. For this we 
note that performing computations on S^xR would be another way to impose an IR cutoff. 
In that situation it would seem that a single gluon would correspond to a configuration 
containing spiky strings, similar in spirit to the ones in ||3^ . It might be possible to consider 
the scattering of spikes by relating them to some objects on the spin chain. 

Hopefully, these amplitudes might be helpful for understanding aspects of QCD and 
the transition between the perturbative and non-perturbative regimes. 
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5. Appendix A: Solution in the original global coordinates 

We start by considering the solution in terms of the original Poincare coordinates (we 
set i? = 1 as it doesn't play any role in the following discussion.) 

a;o = — V 1 + 6^ (cosh^ W2 — cosh^ ui) , 2 = - (coshui cosh W2 -|- fesinhui sinhu2) , 
ZiCi a 

i f U2 1 , ^ 1 , '^1 1 \\ 

xi =- ( ^ + ^ smh2'U2 + o( — ^ + ^ smh2ui) I , 

i f Ul 1 . 1 <^ , /1t2 1 ■ 1 r. \ 

X2 [ smh2'Ui + o[ smh2tt2) 

(5.1) 

In order to study the intersection of our surface with the boundary of AdS we write the 
solution in term of embedding coordinates 

X'^ = — , // = 0,---,3 

X - ^ X X _^!±^ 
A_i -|- A4 — — , A_i — A4 — 

Z Z 

The boundary of AdS in global coordinates is then parametrized by the space of coordinates 
satisfying 

-Xii-X^+Xi2 + ... + X| = (5.3) 

quotientcd by overall rescalings, but they cannot all be zero. The geometry of the boundary 
is i? X S^. 

As discussed in the main text, we impose boundary conditions at z = 00. Notice that 
z = 00 implies the equation X_i -|- X4 = which gives a lightlike surface on the bulk 
that intersects the boundary also on a lightlike surface. The region on the boundary that 
is bounded by this lightlike surface is conformal to Minkowski space. In fact, that region 
corresponds to the Minkowski space parametrized by x'^ at 2; = 0. However, the z = 00 
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surface also extends in the interior of AdS^. In order to see whether the region ±00 is 
on the boundary or not we need to check the behavior of the x'^ coordinates. Without loss 
of generality lets assume mi ^ M2 ~ 1. In this regime we obtain the following expression 
for the embedding coordinates!! 



^ _ iVl+¥ ^ _ ib e"^ ^ 

° 4 cosh U2 +b sinh ^2 ' ^4 cosh U2 + b sinh U2 ' ^ 

^2 — — T ; ~y — r-; : ^-1 — 



4 cosh U2 + b sinh U2 8a (cosh U2 + b sinh W2 ) 

(5.4) 

These coordinates are complex but large. Thus we can say that they live in the complexified 
boundary, which is defined by ( ^.3|) but now the coordinate are complex and we allow 
complex rescalings. 



6. Appendix B: A useful integral. 

In this appendix we show how to compute the following integral 

/"OO 

S= (cosh u cosh + 6 sinh w sinh t;)"^ (6-1) 



This integral can be done by expanding the integrand as a power series on b, then inte- 
grating term by term and finally performing the sum 

^ = f2 r[^+t-\]n ^'' (/ Mcoshunts^nhuyy (6.2) 
The double integral becomes a sum of integrals on a single variable 

duicoshuYitanhuy = + ^".r|!+!l!j^^ — (6-3) 
This identity is valid when e < 0. Finally, performing the sum we obtain 

This expression is valid for all values of e < 0. However, we are interested in the behavior 
for small e. We find that 

2Fi(^,-|,^;62) = l + ^log(l-62)6+llog(l-6)log(l + 6)e2 + 0(e3) (6.5) 
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We thank Shesansu Pal for pointing out an error in the original version of these equations. 
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We find it convenient to express this expansion as follows 



*y 1;"/ (6.6) 
= i + 1 log(l + >>) + j(log(l + l>)f = i(l + b)', F< = i(l - b)' 

Where all the equalities are valid up to order e^. 



7. Appendix C: Extrapolating the weak couping results for g(X) to strong 
coupling. 

In this appendix we give a rough estimate for g{X) at strong coupling by using its 
weak coupling expansionil. We follow closely the idea proposed in in order to obtain 
an approximate expression for the cusp anomalous dimension valid at all values of the 
coupling constant. 



The weak coupling expansion for the cusp anomalous dimension is [39 



A \ / A \^ / A " ^ 



This perturbative information together with its strong coupling behavior /(A) ~ -\/A sug- 
gest a relation of the form 

=T.Cr{f{W (7.2) 



For a given value of n, once we fix the coefficients by using the perturbative data, ( |772 
predicts a given strong coupling limit, il 

At weak coupling g{X) was computed up to three loops 



A V / A 



^W = -4C3(y^) +8(4C5 + 10/3C2C3)(^) +... (7.3) 



This calculation was suggested to us by Z. Bern, L. Dixon and R. Roiban , who also did it. 
Actually, see discussion in p^ , if one introduces enough perturbative information, namely 
n > 4, the value predicted for /(A) at strong coupling is in very good agreement with the known 
strong coupling value. Furthermore, combining the known weak and strong coupling information 
one can get an answer accurate to 1 % (when compared against the answer from the BES equation) 
for the whole range of the coupling constant. 
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To repeat the analysis discussed above in order to extract the strong coupling behavior 



of g{X) is subtle due to two features. First, g{X) starts at two loops, so the relation ( |77^ ) 
should be modified accordingly. Second, note that if our prediction is correct, g{X) should 
change its sign in going from weak to strong coupling. This does not occur for /. 

A possibility is to define a new function ^(A), differing from g{X) by a constant times 
/(A). As seen in the main text, this shift corresponds simply to a change in the IR 
regulator. 



~g{X)=giX)+^fiX) (7.4) 

For positive ^, ^(A) will be positive both at weak and strong coupling. It is then 
straightforward to repeat the analysis of for n = 3 (since we have three coefficients at 
our disposal). The value obtained at strong coupling should be then compared with our 
prediction 



^(A) = (2(l-log(2))+40^ + ... (7.5) 

The following figure shows a comparison between the value predicted by our compu- 
tation ( [f75|) (solid blue line) , and the value obtained by a naive extrapolation as explained 
above (dashed red line), for the range < ^ < 2 




Fig. 6: Comparison between the result for g obtained from our computation (solid 
blue line) and a naive extrapolation from weak coupling (dashed red line). 



Note that for ^ = the assumption (|7.2| ) is not valid. At the particular value ^ = 
log(2)/2 the strong coupling prediction is particularly simpleHl 

~g{X) = 2^ + ... (7.6) 
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Note that if log(2) has transcendentality 1, then such a shift doesn't break transcendentality. 
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whereas the value from the extrapolation turns out to be 

^(A)^1.37^ + ... (7.7) 

This differs by a 30% from the predicted number. All together, given the little that is 
known for the function g{X), all that we can say is that a naive extrapolation from weak 
to strong coupling seems to be compatible with our results. 

Of course, the above analysis is by no means rigorous. In order to make a more precise 
comparison one would like to have the analog of the BES equation from which g{X) could 



be computed at any order, then one could proceed along the lines of [40,41 
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